Magnetic Properties of (VO)2P207 from Frustrated Interchain Coupling 
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Neutron-scattering experiments on (VO)2P2C>7 reveal both a gapped magnon dispersion and an un- 
expected, low-lying second mode. The proximity and intensity of these modes suggest a frustrated 
coupling between the alternating spin chains. We deduce the minimal model containing such a 
frustration, and show that it gives an excellent account of the magnon dispersion, static suscepti- 
bility and electron spin resonance absorption. We consider two-magnon states which bind due to 
frustration, and demonstrate that these may provide a consistent explanation for the second mode. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.40. Gb 

Vanadyl pyrophosphate ((VO)2P2C>7, abbreviated 
VOPO)u is a low-dimensional quantum magnetic system 
composed of S = \ V 4+ ions. These have Heisenberg an- 
tifenjomagnetic (AF) interactions, and a singlet ground 
stated with spin gap A = 3.1meV. Based on suscepti- 
bility and neutron-scattering experiments on powders, 
VOPO has been compared in detail to spin-ladder and 
alternating-chain models ,□ both of which were found to 
be consistent with the data. Rising interest in spin lad- 
ders, along with the assumption that the strongest ex- 
change paths would be those between V5+ ions with the 
smallest separations, led to a preferenceM for the ladder 
conformation, and for several years VOPO was quoted as 
the first known spin ladder system. 

This picture was conclusively debunked by Garrett and 
coworkers, who performed the first neutron-scattering 
measurements on aligned single crystals of VOPO.tl 
These demonstrated unequivocally that the strongest ex- 
change path in the system ( J x ) was a double V-O-P-O-V 
link through phosphate groups in the crystallographic in- 
direction, while the next (J2) was a double V-O-V link 
between edge-sharing VO5 square pyramids, also along 
b. The V-O-V bond along a through the apex of the 
pyramid, believed to be the strong ladder leg, was found 
to be very weak. The strength of the V-O-P-O-V bond 
was verified independently! 4 ] in the related compound 
VODPO4. iD 2 0, where it appears in isolation. These re- 
sults are not unexpected, because the single electron on 
V 4+ in a square pyramidal environment occupies the d xy 
orbital in the basal (be) plane of the unit. They lead to 
the interpretation of VOPO as a set £>f dimerized chains, 
with coupling ratio A = J2/J1 — 0.8fl 

The same experiment obtained the most detailed mea- 
surements to date for a second mode withj-gap 5.7meV. 
This was identified from a previous studyO as a candi- 
date two-magnon bound state. In addition, the inter- 
chain coupling was deduced to be weakly ferromagnetic 
(FM), a result unexpected for a conventional, lwg su- 
perexchange path. We begin from the observationErLj that 
frustration can act to promote bound states between ex- 
citations, and the hypothesis that the FM interchain cou- 
pling in fact results from a competition between two AF 



couplings. We will deduce and justify the nature of such 
a frustrated, two-dimensional (2d) model for VOPO, and 
show that it provides the basis for a complete account of 
the magnetic properties of this material. 

The structure of VOPO was determined in detail by 
Nguyen et alia The 6-axis alternating chains are coupled 
along a by either the V-O-V path through apical O, or 
by V-O-P-O-P-O-V pathways through two PO4 groups, 
while c-axis coupling (see also Ref. |l|) proceeds through a 
single phosphate unit. Quantitative superexchange cal- 
culations, particularly for extended pathways, remain be- 
yond the scope of current understanding and computer 
power particularly for the V ion. We adopt here a qual- 
itative approach of identifying paths coupling strongly to 
the V d xy orbital, and using consistency with experiment 
to construct the minimal set of necessary interactions. 
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FIG. 1. Schematic representation of the VOPO system, 
showing frustrated coupling between dimerized chains. 

We discard the possibility that the V-O-V bond along 
d is significant, because its direction is orthogonal to the 
d xy orbital. We argue instead in favor of the V-O-P-O-P- 
O-V path, due to the demonstrated importance of phos- 
phate groups, and because bond angles throughout its 
length can remain close to 180°. Such paths connecting 
V ions separated only along d (J a ), and those separated 
both along a and by one ionic spacing along the chains 
(Jk), are in principle rather similar, which presents the 
possibility of frustration. However, there are 2 types of 
V-O-P-O-P-O-V path, those connecting into the phos- 
phate groups mediating J\ and those connecting to the 
O atom concerned in J2. If the latter were significant, 
one could expect in addition a strong c-axis coupling, 
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which is not observed.i We assume therefore that path- 
ways of the second type are strongly suppressed by the 
angles and orbital configurations in the V-O-P overlap, 
and concentrate on the resulting 2d model in Fig. 1. 

Here J\ = J is set to 1, A is defined above, and 
(J-a = Ja/J and fib = Jb/J parameterize the competing 
V-O-P-O-P-O-Vsuperexchange interactions. We exclude 
a second-neighbor intrachain coupling, due to both the 
lack of a suitable exchange pathway and the minimacy 
criterion, as it is found not to be important to any of the 
quantities we will compute. That fi a ib maypbe significant 
is shown by a high-temperature expansionH of the static 
susceptibility x(T)- This gives the Curie temperature 9 C 
as the sum of the interactions of a structural unit accord- 
ing to 9 C = + L+2n a + 2fib]. Comparison with the 
measured 6 C = — 84KL!l suggests that the chain couplings^ 
alone satisfy at most 75% of the sum. 

We analyze this model using the method of Ref. ||. 
Taking the strongest (Ji) bonds as the sites of dimers, 
whose singlet ground states may be excited to triplets, 
the interactions permit delocalization, or hopping, of the 
triplets. The hopping gives a kinetic energy which lowers 
the bond singlet-triplet gap to the physical one, and a 
description of a dispersive, triply-degenerate 1-magnon 
excitation in 2d reciprocal space. While this technique is 
best suited to strong]*' dimerized systems, perturbation- 
theoretic calculationsB show rapid suppression of higher- 
order contributions even for rather weak dimerization, 
and thus internal consistency. 

We have calculated the coefficients for hopping of a 
triplet excitation by perturbative expansion to 5th order 
in the parameters of Fig. 1. Contributions beyond 3rd 
order indeed remain small (~3%), and we illustrate the 
method to this order. Let be the coefficient for the 
excited triplet to hop by i dimer bonds along the chain, 
and across j chains, and fi± denote fi a ± fib, then 



(Ji, J 2 , J a , J b ) = (10.7, 8.49, 2.17, 2.73)meV. 
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The mode dispersion is given most straightforwardly by 

w(q) = J ^2 2V- 5 ™- s i°h ij cos(iq y ) cos(jq x ), (2) 
ij 

although in fact to this order one obtains a better fit by 
expanding the (smoother) quantity w 2 (q). 

The measured dispersion data for the lowest-lying 
mode in Ref. combined with the condition on 8 C , al- 
low one to fit the optimal model parameters as the set 
(J, A, (i a , jib) = (10.7meV, 0.793, 0.203, 0.255), which cor- 
respond to the superexchange interactions 



The various fits offer minor differences in the exact shape 
of the dispersion curve, with the higher-order calculations 
returning a very good description of the overall form. For 
clarity we show only the 5th-order fits in Figs. 2 (a) and 
(b) for the reciprocal-space cuts where data is available. 
The results (§) may be used to predict the form of the 
1-magnon dispersion where it has not yet been measured, 
and this is shown in Figs. 2(c) and (d). We expect the 
upper band edge to be observed at 15.9meV. 
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FIG. 2. Magnon dispersion (solid lines), continuum edge 
(dashed lines) and bound S = 1 states (dotted lines) for (a) 
(0,2tt - q), (b) (q,2n), (c) (vr, 2vr - g), and (d) (2tt - q,n). 
Symbols are data taken from Ref. [3]. 

The primary features of the optimal parameter set are 
as follows. The alternating chains have the rather weak 
dimerization A ~ 0.8, or S = (1 - A)/(l + A) ~ 0.11. The 
interchain couplings are significant, with a magnitude 20- 
25% of the primary chain energy scale. The cross cou- 
pling Jb is slightly larger than the pure a-axis coupling 
J a , which is required for a net FM interchain dispersion, 
and fully plausible from above Fig. 3 shows the density 
of states determined from the dispersion. The dominant 
feature is the logarithmic singularity at the saddle point 
(tt, 0), with energy 4.55meV. This energy, and not simply 
the minimum gap A = 3.1meV at (0,0), will play a sig- 
nificant role in determining thermodynamic quantities. 

We next consider two further magnetic properties for 
which independent measurements are available, to deter- 
mine the consistency of the model and parameters. The 
static, uniform susceptibility x(T) was first measured fat 
powder samplcs,EI and more recently for single crystals.EHI 
The susceptibility to applied field H measures available 
excitations with AS = 1 , so is primarily a probe of the 1- 
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magnon branch. It is given in general by x = —d 2 F/dH 2 , 
where F = — \nZ is the free energy and /3 _1 = ksT . 
In the dimer model, an excited triplet is effectively a 
hard-core boson, because only one such state is physi- 
cally possible on each bond. Thus thermodynamic cal- 
culations require a correction from Bose statistics to ex- 
clusion statistics appropriate for triply degenerate exci- 
tations. The free energy per dimer isliil 



/ = -/T 1 ln{l + [1 + 2 cosh(/3h)]z(P)}, 



(4) 



where h denotes g^sFl and z((3) is the partition function 
for a single mode. The susceptibility per site is then 



Xo= MS)/ (1 + 3*08)): 



(5) 



in which the denominator suppresses the Bose result at 
high T by excluding multiple mode occupation. 
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FIG. 3. Density of states for computed magnon dispersion. 

While the statistical factor corrects for state availabil- 
ity, it does not take account of magnon interaction effects. 
These we may include within the following mean-field 
scheme. Approximating any interaction term with a spin 
on a neighboring dimer by J^+sSj-Sj+a — > Js(S-)S^ +s = 
J$mSfi, the instantaneous magnetization per site m con- 
tributes to an effective internal magnetic field 




hint = —mJA — 2mJ(fi a + /Z&) = —mC. 



(6) 



The susceptibility is defined by m = x^ext, while the 
magnetization m = Xohtot, where h tot = h ext + h mt is 
the total field at each site. Simple rearrangement yields 



X = xo/(l + C X o), C = J(A + 2 M+ ), 



(7) 



as the mean-field, interaction-corrected susceptibility. 

In Fig. 4(a) is shown the quantity x(T) (@). The 
agreement with both powder and single-crystal data is 
good. Qualitatively, both the exclusion statistics factor 
and (Fig. 4(b)) the mean-field correction are required 
to deduce x(T) within this framework. Quantitatively, 
the temperature T max of the peak in the model is 58K, 
rather lower than both data sets, while Xmax is fraction- 
ally smaller. This latter result corroborates the presence 
of frustrating interactions, whose mean-field correction 
gives a suppression. From Fig. 4(b), the model 6 C is in 
excellent accord with the data. 




400 500 

T[K] 

FIG. 4. (a) x(T\-£fom model parameter set. Deviations 
from measured datal_l are discussed in text, (b) X _1 (^) with 
(solid line) and without (dashed) mean-field correction. 

The discrepancies between calculation and experiment 
may be analyzed by applying the triplet hopping theory 
with mean-field correction to the dimerized chain. For a 
chain with A = 0.8, the resulting x(T) is directly cg 
parable with essentially exact numerical simulationsJ13 
This comparison (with Figs. 2 of Refs. |lj and ||) reveals 
that T^ x /T^ ax ~ 0.75, while xZaJx n max - 1-2, where 
superscripts m and n denote model and numerical results 
respectively. The first ratio is very similar to the discrep- 
ancy between model and experiment above. The second 
suggests further that in Id the mean-field correction is 
not sufficient to reproduce interaction-induced suppres- 
sion of x, whereas improved agreement is achieved in 2d. 
We conclude that deviations from the data arise primar- 
ily because the calculational scheme cannot fully account 
for magnon interactions, and not due to any intrinsic 
shortcomings of the model. 

ESR experiments have also been performed on pow- 
der and crystalline VOPO samples, and we focus on 
the latter .E3 For linearly-polarized (microwave) radiation, 
which excites-transitions of A5 Z = ±1, the power absorp- 
tion is givenlij by the imaginary part of the susceptibility 



x"(w) = nS(huj — h) smh((3Tiw) 



d 2 k 



-/Mk) 



(8) 



Including the hard-core boson constraint as above, the 
intensity at the resonance frequency Tilo = h is 



7(/3)oc 



sinh.(J3h)z(0) 



l + (l + 2cosh(/%))z(/3)' 



(9) 



This quantity contains in principle the mean-field cor- 
rection in the form h = h ext — Cm. However, for the 
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strongest resonance at 134GHz ~ 7K, the sinh function 
is linear at temperatures on the order of A, and alter- 
ations of h affect only the prefactor. 
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FIG. 5. ESR absorption from model parameter set, show- 
ing excellent consistency with observation £j 



The calculated ESR absorption at 134GHz is shown in 
Fig. 5 with the data of Ref. [l0[ The excellent agreement 
indicates the importance of the 2d density of states, in 
particular the energy of the logarithmic singularity (Fig. 
3) . While the ESR absorption is often considered to mea- 
sure a q = property akin to x(T), we note here a signifi- 
cant difference between the two in experiment. Theoreti- 
cally, the triplet hopping approach yields a single energy 
scale for both quantities, which is in more satisfactory 
agreement with the ESR absorption, indicating that this 
is rather less sensitive than xCO to magnon interactions. 

We conclude with a brief analysis of possible bound 
states of two magnons with the computed dispersion, 
to seek an explanation for the observed second mode. 
The initial perturbative problem, represented generally 
as H = Hq + aH' , may be brought by a continuous, uni- 
tary transformationt-3 to a new Hamiltonian H c s which 
conserves triplet number ([H e s , Hq] — 0). This procedure 
is conducted perturbatively in a to generate a series for 
-ffcff- The action of H e s on the sector with one triplet 
gives the one-magnon dispersion, while on the two-triplet 
sector with total spin S = 1 it yields the two-magnon con- 
tinuum and contains in addition magnon interactions. At 
fixed momentum one obtains a 2-body problem soluble 
by Lanczos tridiagonalization, in particular for the en- 
ergy of (triplet) bound states. The number of interac- 
tion coefficients rises very rapidly with order n in a. We 
have completed calculations for n — 3 and 4, finding only 
small alterations (« 2% of the continuum edge energy), 
but caution that these remain potentially significant. 

The predicted dispersion curves for the two-magnon 
bound state, computed at 4th order, are shown as the 
dotted lines in Fig. 2. We find that with the parame- 
ters fixed as above, the bound state does appear as an 
excitation mode below the continuum over a large part 
of the Brillouin zone. However, this is not the case close 
to the band minimum, where it is lost in the continuum. 
While we thus do not always find a second mode below 
2A in a consistent treatment, the discrepancy between 
the computed bound state and the observed resonance 



is not large: it is reasonable to postulate that forthcom- 
ing refinement of this calculation, and possible extension 
of the model to include further interactions, will indeed 
yield a quantitative explanation of the second mode. 

In summary, we have presented a model for frustrated, 
2d coupling in (VO)2P20y. The model gives an excellent 
and fully consistent account of the available data concern- 
ing elementary excitations, namely the 1-magnon disper- 
sion, static susceptibility and ESR absorption. We use it 
to predict the dispersion in the entire Brillouin zone, and 
to indicate the origin of the observed, low-lying second 
excitation mode as a triplet 2-magnon bound state. 

During completion of this woek, we became aware of 
the contribution of Weisse et aZ.E-3 These authors investi- 
gated the same model, albeit without the physical justi- 
fication presented here, by diagonalization on small (up 
to 4x8) clusters. Their results are qualitatively in good 
agreement with our expectation that frustration should 
promote bound states. Quantitatively, the parameters 
chosen differ considerably, and we would not expect to 
see good agreement with susceptibility and ESR data. 

We are grateful to U. Low, B. Liithi, S. Nagler, H. 
Schwenk, and particularly D. A. Tennant for invaluable 
discussions and provision of data. We thank also E. 
Muller-Hartmann for helpful conversations, and C. Knet- 
ter for assistance with the calculation of H c s. GSU was 
supported by an individual grant and by SFB 341 of the 
DFG. BN wishes to acknowledge the generosity of the 
Treubelfonds. 



1 D. C. Johnston, J. W. Johnson, D. P. Goshorn, and A. J. 
Jacobson, Phys. Rev. B 35, 219 (1987). 

2 T. Barnes and J. Riera, Phys. Rev. B 50, 6817 (1994). 

3 A. W. Garrett, S. E. Nagler, D. A. Tennant, B. C. Sales, 
and T. Barnes, Phys. Rev. Lett. 79, 745 (1997). 

4 D. A. Tennant, S. E. Nagler, A. W. Garrett, T. Barnes, 
and C. C. Torardi, Phys. Rev. Lett. 78, 4998 (1997). 

° T. Barnes, J. Riera, and D. A. Tennant, preprint # cond 



mat 9801224 Phys. Rev. B (to appear). 
G. S. Uhrig and H. J. Schulz, Phys. Rev. B 54, R9624 

(1996); Phys. Rev. B 58, 2900 (1998). 
7 G. Bouzerar, A. P. Kampf, and G. Japaridze, preprint # 



cond- mat 9801046, Phys. Rev. B 58, 3117 (1998). 



8 P. T. Nguyen, R. D. Hoffman and A. W. Sleight, Mat. Res. 
Bull. 30, 1055 (1995). 

9 G. S. Uhrig, Phys. Rev. Lett 79, 163 (1997). 

10 A. V. Prokofiev, F. Biillesfeld, W. Assmus, H. Schwenk, 
D. Wichert, U. Low, and B. Liithi, Eur. Phys. J. B 5, 313 
(1998). 

11 M. Troyer, H. Tsunetsugu, and D. Wiirtz, Phys. Rev. B 
50, 13515 (1994). 

12 J. C. Bonner, H. W. J. Blote, J. W. Bray, and I. S. Jacobs, 
J. Appl. Phys. 50, 1810 (1979). 

13 I. Affleck, Phys. Rev. B 46, 9002 (1992). 



4 



14 p 

15 A. 



J. Wegner, Ann. Physik 3, 77 (1994). 

Weisse, G . Bouzerar, and H. Fehske, preprint # cond- 



mat 9805374 



5 



